Apparent Thermalization due to Plasma Instabilities in Quark-Gluon Plasma 
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Hydrodynamical modeling of heavy ion collisions at RHIC suggests that the quark-gluon plasma 
(QGP) "thermalizes" in a remarkably short time scale, about 0.6 fm/c. We argue that this should 
be viewed as indicating fast isotropization, but not necessarily complete thermalization, of the non- 
equilibrium QGP. Non-Abelian plasma instabilities can drive local isotropization of an anisotropic 
QGP on a time scale which is faster than ordinary perturbative scattering processes. As a result, we 
argue that theoretical expectations based on weak coupling analysis are not necessarily in conflict 
with hydrodynamic modeling of the early part of RHIC collisions, provided one recognizes the key 
role of non-Abelian plasma instabilities. 



Hydrodynamic models of RHIC collisions (based on 
near-ideal fluids) provide a good description of a wide 
range of experimental data, including radial and ellip- 
tic flow measurements, provided one assumes that the 
initial partons thermalize in about 0.6 fm/c [1]. How- 
ever, theoretical estimates based on perturbative scat- 
tering processes yield expected thermalization times in 
the range of 2.5 fm/c or above [2]. What is the signifi- 
cance of this discrepancy? Are weak-coupling analyses, 
which should be valid for asymptotically high energy den- 
sities (and asymptotically large nuclei), inapplicable at 
RHIC energies? Perhaps so. Or have dynamical pro- 
cesses which may be responsible for this fast apparent 
thermalization not been correctly identified? We will ar- 
gue this is the case [3]. Estimates based on perturbative 
scattering neglect essential dynamics: the collective be- 
havior associated with non-Abelian plasma instabilities. 
Such instabilities can produce large non-perturbative ef- 
fects, including apparent thermalization. We discuss two 
qualitative lessons which emerge from a weak-coupling 
analysis: (i) hydrodynamic behavior does not require full 
thermalization — isotropization of parton momenta in lo- 
cal fluid rest frames suffices; and (ii) plasma instabilities 
can drive isotropization at rates which are parametrically 
faster than perturbative scattering rates. 

APPARENT THERMALIZATION 

The thermalization time scale in a quark-gluon plasma, 
defined as the inverse relaxation rate of arbitrarily small 
departures from equilibrium, depends on the rate of 
large-angle scattering (and near-collinear splitting/join- 
ing) processes among quarks and gluons [4]. Parametri- 
cally, this time scale is [5] ~ 26 [g A T\n(2A/ g)]- 1 , and for 
plausible values of RHIC parameters it is hard to recon- 
cile this time scale with the fast apparent thermalization 



observed in RHIC collisions. However, this time scale 
characterizing relaxation of asymptotically small pertur- 
bations is irrelevant to the question of when hydrody- 
namic models can be a good approximation to the dy- 
namics of a non-equilibrium quark-gluon plasma. The 
essential assumption of ideal fluid hydrodynamic models 
is that the stress tensor, in the local rest frame at some 
point in the system, is nearly diagonal, 
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with some equation of state relating the pressure p to 
the energy density. But relation (1) is just a statement 
of isotropy (in the local fluid rest frame), and is automat- 
ically true if typical excitations have random directions 
— even if their energy distribution is far from thermal, 
or if the pressure p differs from the equilibrium pressure 
for a given energy density. Consequently, understanding 
when a hydrodynamic model can first provide a good ap- 
proximation to the plasma dynamics is the same question 
as understanding what dynamics drives isotropization. 

PLASMA INSTABILITIES 

To begin, we summarize known results concerning 
gauge field instabilities in anisotropic non-Abelian plas- 
mas. Further details may be found in Refs. [4, 6-8]. 

Let Phard denote the characteristic momenta of typi- 
cal excitations in a non-equilibrium quark-gluon plasma. 
(For example, in the saturation scenario [9], phard equals 
the saturation scale Q s at time Qj 1 -) We assume that 
Phard is sufficiently large that these excitations act like 
highly relativistic particles. For time scales short com- 
pared to the mean free time between large-angle scatter- 
ings of typical excitations (and large compared to p^ard)' 
the natural framework for describing the dynamics is col- 
lisionless kinetic theory. One splits the degrees of freedom 
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into short wavelength (or "hard" momentum) excitations 
which may be characterized by a phase space distribution 
function /(p, x, t), and long wavelength (or "soft" ) gauge 
field modes which may be regarded as forming a classical 
field. For a non-Abclian theory, the resulting Boltzmann- 
Vlasov equation has the form [10, 11] 

(A + v • Ac) / + \9 {(E + v x B) it V Pi /} - . (2) 

The corresponding Maxwell equations are 



(D,n a =j:^/ P ""tr(u): 
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with ,f p = / ( 2 7r)9 ' uA1 = (1'P)' an( i *o a c °l° r generator. 

Any distribution which is homogeneous (in space) and 
colorless, combined with vanishing soft gauge field, gives 
a static solution to Eqs. (2,3). Perturbations about such 
solutions obey a linearized equation of motion [obtained 
by linearizing Eq. (2) in deviations from the static solu- 
tion, solving for Sf, and plugging the result into Eq. (3)] 
which (after a space-time Fourier transform) has the form 



{K 2 g» v - K»K V + W V {K)} A V {K) =0. 



(4) 



where the wavevector K 1 * = (w,k) [12]. The retarded 
gauge-field self-energy, generated by hard excitations, is 



W v {K)=g 2 [ 



am 

p dp 1 



v»g lv + 



v ■ K — ie 



(5) 



The zero-frequency spatial self-energy 11^(0, k) depends 
on the direction but not the magnitude of the spatial 
wavevector k. If /(p) is anisotropic but parity invariant, 
then the self-energy matrix 11(0, k) has a negative eigen- 
value for some directions of k. This implies that there 
are unstable solutions to the small fluctuation equation 
(4), i.e., solutions for which oj has a positive imaginary 
part [4, 6]. These are non-Abelian versions of Weibel 
instabilities in ordinary plasma physics [13]. 

Let — /i 2 denote the most negative eigenvalue of 11(0, k) 
(for any k). Unstable modes have |k| < \i. Let 7 de- 
note the maximal growth rate of unstable modes. If the 
hard particle distribution has 0(1) anisotropy [14] then 
the maximum unstable wavevector fj, and the maximum 
growth rate 7 are both comparable to the effective mass 
TOqo of hard gluons, 
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If Phard is the momentum scale which dominates the in- 
tegral (6), and n = J p /(p) is the spatial density of hard 

excitations, then m.00 ~ g y ^/Phard- 

To compare to perturbative scattering rates consider, 
for example, a system with n = 0(pjj ard ) — the same 
parametric relation as in equilibrium, where p ~ T 



and n = 0(T 3 ). In this case TOto, and hence the in- 
stability growth rate 7 for 0(1) (or larger) anisotropy, 
is 0(g phard)- This rate is parametrically faster than 
the 0(g 4 Phard) rates for large-angle scattering or near- 
collincar splitting, or even the 0(g 2 Phard) rate of small- 
angle scattering [4]. More generally, for 0(1) anisotropy 
7 is faster than the large-angle scattering rate whenever 
n <C Phard/ff 2 [15]- This inequality is satisfied paramet- 
rically unless there is saturation, and even in saturation 
scenarios, it is satisfied for t ^> Qj 1 [9]. 

Numerical values depend, of course, on the specific 
form of the anisotropic phase space distribution. A sim- 
ple example [16] involving a typical particle energy of 
1 GeV, plasma energy density of 27 GeV/fm 3 , a phase 
space distribution proportional to (p • z) 4 , and a s = 0.5 
yields ~ 740 MeV, and 7 ~ 280 MeV = (0.7 fm/c)" 1 
for k ~ 575 MeV. With more extreme anisotropy, the 
growth rate 7 can approach moo itself [6] . Yet other an- 
gular distributions can give slower growth rates. 

Instabilities will grow exponentially until some dynam- 
ics comes into play which causes the amplitudes of unsta- 
ble modes to saturate. There are two natural possibilities 
for when this might happen [17]. If the unstable modes 
with wavenumbers of order fi grow until the soft gauge 
field has an 0(n/g) amplitude [or the field strength is 
0(fj, 2 /g)}, then non-Abelian corrections to the linearized 
equation of motion (4) will become important and could 
substantially affect the further evolution [18]. In partic- 
ular, one might expect these non-linearities to lead to 
efficient transfer of energy from the unstable modes to 
stable modes (with comparable wavenumber). 

Alternatively, if instabilities do not saturate at 0(fj,/g) 
amplitudes, then they may continue growing until their 
amplitudes reach the scale phard/ff [and field strengths 
are 0(/xphard /<?)]■ This is the point where the soft gauge 
field no longer acts as a small perturbation on the mo- 
tion of hard excitations. To see this, note that for this 
amplitude, the gauge field part of a covariant derivative 
is just as large as the ordinary derivative when acting 
on fluctuations with O(phard) momenta. This is also the 
point where the energy density in the soft gauge field 
becomes an 0(1) fraction of the total energy density, 
( F soft) 2 ~ (^Phard/5) 2 ~ np hald . 

There are reasons to believe the second alternative, not 
the first, is correct. The generalization to anisotropic 
plasmas of the HTL ("hard thermal loop") effective ac- 
tion is [19, 20] 
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Evaluating this, explicitly, for arbitrary static fields in 
order to examine the corresponding effective potential is 
not feasible. But in the special case of fields which vary 
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in only one spatial direction, the effective action reduces 
to a simple local form. Let n denote the direction of the 
wavevector of the most unstable mode. For gauge fields 
which depend only on n • x, one finds that the effective 
potential is [21, 22] 

V[A(n-x)] = J Sx [\F9.F^ + $A$H ij (0,&)Afi. (8) 

When 11(0, n) has a negative eigenvalue this potential is 
unbounded below. The runaway directions of steepest- 
descent correspond to Abelian field configurations where 
the commutator terms in the field strength vanish. 
This suggests that non-Abelian non-linearities may not 
cause growing instabilities to saturate at the scale n/g, 
provided the field configuration evolves toward an effec- 
tively Abelian form which can continue rolling down the 
potential energy landscape. This behavior has been seen 
in time-dependent numerical simulations in 1+1 dimen- 
sions [22, 23] — the instability locally "Abelianizes" and 
continues growing. It is important to perform full 3+1 
dimensional simulations of the collisionlcss kinetic the- 
ory (2,3) to verify this conclusion. Such simulations are 
in progress [24]. Here, we shall assume that growth of in- 
stabilities, beyond the soft scale fi/g, will be confirmed. 

ISOTROPIZATION 

Growing instabilities imply that the stress tensor of 
the non-equilibrium system will receive growing contri- 
butions from the soft gauge field. The fastest growing lin- 
earized modes tend to decrease the anisotropy in the to- 
tal stress tensor [7]. For example, if the anisotropic hard 
particle distribution has a prolate form, so that T^f rd 3> 
T^, T^ rd , then the wavevectors of the fastest growing 
unstable modes lie in the equatorial plane and the growth 
of these modes produces a soft gauge field contribution to 
the stress tensor which is oblate, TJ° ft ~ T y s ° ft > T|° ft . 
Conversely, for an oblate hard particle distribution, the 
fastest growing unstable mode has its wavevector along 
the normal direction and generates a prolate contribution 
to the stress. Hence, even in the linearized regime, one 
can see that soft gauge field instabilities push the system 
toward greater isotropy. However, the soft contribution 
to the stress tensor is small compared to the hard particle 
contribution, and the back-reaction of the soft gauge field 
on the hard particles is a tiny perturbation, as long as the 
soft gauge field amplitude is much less than O(phard/s0- 

But if the soft gauge field amplitude reaches the scale 
Phard/<7 then it no longer acts as a small perturbation 
to the dynamics of hard excitations. Recall that the ra- 
dius of curvature of an excitation of momentum p and 
charge g in a magnetic field B is R = p/(gB). If the 
radius of curvature is comparable to the magnetic field 
coherence length which means B ~ [i-p/g, then ex- 
citations of momentum p will undergo O(l) changes in 



direction during traversals of any single coherence-length 
sized magnetic field "patch" [25]. 

Therefore, if unstable soft gauge field modes with 0(/z) 
wavevectors grow until the field strength is 0(fiphard/g), 
then typical excitations will experience 0(1) changes in 
direction in times of order pT 1 . Excitations with differing 
momenta or colors will receive different deflections from a 
given patch of (non-Abelian) magnetic field. Excitations 
traversing different patches of magnetic field (separated 
by 0(/x -1 )) will receive nearly uncorrelated deflections. 

The net effect is that a soft gauge field with a non- 
perturbative amplitude of order Phard/g can effectively 
drive isotropization in the distribution of typical hard 
excitations on a time scale which equals the coherence 
length of the soft gauge field. And isotropization of 
the hard particle distribution will turn off further growth 
in the soft gauge field (since gauge field instabilities are 
absent for isotropic distributions). 

As with all instabilities, the time, or number of e- 
foldings, required for the soft gauge field to become large 
depends on the size of initial "seed" amplitudes in the rel- 
evant unstable modes. The amplitude of the soft (k ~ /i) 
gauge field generated by a random color charge distri- 
bution of the hard particles can be estimated as A 2 <~ 
g 2 n//i ~ g^nphmd- This is the smallest the seed field 
could be. For densities from n = 0(pf 1&!:d ) U P to the 
density limit n = 0(p 3 lald /g 2 ) imposed by saturation, 
A ^ 0(g 1 ' 2 phard)- This is only a factor of g 3 ' 2 smaller 
than the nonperturbative 0(phard/ <?) amplitude. There- 
fore, the number of e-foldings required for instabilities to 
grow to this non-perturbative size is only of order ln(l/g) . 

Treating logs of g, for simplicity, as 0(1), this means 
that if the initial anisotropy is O(l) then the charac- 
teristic growth time needed for unstable modes of the 
soft gauge field to reach the non-perturbative amplitude 
Phard/g is only of order 7 _1 . The resulting soft gauge 
field then drives isotropization of the hard particle dis- 
tribution on a comparable 7 _1 time scale. Therefore, 
(up to logs of g and factors of order one), the time scale 
for isotropization of the hard particle distribution is the 
same as the (inverse) instability growth rate 7 _1 [26]. 

In numerical simulations of ordinary non-relativistic 
plasmas, essentially the same process of instability- 
driven isotropization has been observed [27], with the 
growth of magnetic instabilities driving large reductions 
in anisotropy once the magnetic fields reach critical 
strength. (These simulations allowed three-dimensional 
momentum space variations, but assumed translation in- 
variance in one spatial direction.) Various QGP numer- 
ical simulations [28] have failed to see any sign of this 
instability-driven dynamics because they did not allow 
full three-dimensional variations. 

Although we have focused on the ability of non- 
perturbative soft gauge fields to generate large changes in 
directions of hard excitations, it should be noted that /i -1 
is also the characteristic time scale for 0(1) changes in 
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energies of hard excitations. This is inevitable, given the 
fact, noted earlier, that when the soft gauge field reaches 
the non-perturbative amplitude Phard/ff its energy den- 
sity is comparable to the energy density in the hard ex- 
citations. But it may also be seen directly by noting 
that chromoelectric fields generated during the growth 
of instabilities will be comparable in size to chromomag- 
netic fields (since the growth rate of unstable modes is 
comparable to their wavenumbers for 0(1) anisotropy). 
So chromoelectric fields will reach the same 0(/iphard/<?) 
size as magnetic fields — which means that an excitation 
traveling a distance p _1 will have work of order phard 
done on it by the soft gauge field. Of course, this time 
scale for O(l) changes in energy may be very different 
(and much shorter) than the time scale for true thermal- 
ization, as defined by a near-thermal energy distribution 
of excitations over a parametrically large dynamic range. 

CONCLUSIONS 

We have argued that "early thermalization" in heavy- 
ion collisions is more properly interpreted as evidence of 
fast isotropization in the distribution of excitations. And 
we have argued that non-Abelian plasma instabilities can 
drive isotropization at a rate which is parametrically fast 
compared to perturbative scattering rates. Consequently, 
we see no reason to view the fast onset of hydrodynamic 
behavior in RHIC collisions as necessarily in conflict with 
theoretical expectations based on weak-coupling analysis 
of a quark-gluon plasma, provided one properly accounts 
for the effects of non-perturbative plasma instabilities. 
Further study of the scenario we have sketched is cer- 
tainly needed; in particular full three-dimensional non- 
Abelian Boltzmann-Vlasov simulations with appropriate 
initial conditions should be conducted. 
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